Reverse Engineering with Quantum Noise 



Muhammed Yonag and Joseph H. EberljQ 
Rochester Theory Center, and Department of Physics and Astronomy, 
University of Rochester, New York 14627, USA 

We show that specific quantum noise, acting as an open-system reservoir for non-locally entangled 
atoms, can serve to preserve rather than degrade joint coherence. This creates a new type of long- 
time control over hiding and recovery of quantum entanglement. 
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I. INTRODUCTION 

Usually the "open" part of a universe is associated with 
a reservoir or a noise source and is the part that is poorly 
known or resistant to control. The quantum vacuum is 
the most fundamental example. Here we examine a sit- 
uation in which the opposite is the case, in which the 
open part of the relevant universe, the noise reservoir, 
is known to be quantum radiation and its interaction is 
used to control and manage, actually store safely, pre- 
existing coherence. This is what could be called reverse 
engineering with quantum noise. A case of particular in- 
terest arises when the coherence is shared between parts 
of the system. For example, it has been shown recently 
l] that open quantum system theory has unanticipated 
decoherence features such as non-additive response to 
weak Markovian noise when the coherence is a two-party 
entanglement. Here we demonstrate an example of re- 
verse engineering using quantum noise, the consequence 
of which is a long-term on/off quantum switch for entan- 
glement. 

The effect of radiation fields on entanglement of CQED 
atoms has a long history, beginning with the atom-field 
entanglement exhibited in the original Jaynes-Cummings 
(JC) paper Q. Much later work by Gea-Banacloche Q 
and Phoenix and Knight Q extended these considera- 
tions into the long-time domain that was revealed by the 
discovery of collapse and revival physics d, H, 0]- Still 
later, the evolution of CQED entanglement was further 
extended to the case where the atom and its field start 
from mixed states [8j. 

We note that an entirely new set of opportunities arises 
when one asks about the entanglement of two atoms in 
different CQED cavities. Then the atoms cannot inter- 
act in any way, and pure non-local entanglement enters 
the picture. Examinations of this situation have taken 
into account the back action of the radiation emitted by 
the atoms themselves, and have included treatments of 
the cavities ranging between two extremes: broadband 
absorptive 0, [l(| or mirror- like [ll[ . Among the conse- 
quences was so-called entanglement sudden death (ESD 
- see [1, [H, 0; 0] ) ' which in the mirror-cavity case and 
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its analogs was followed by rebirth, with death and re- 
birth repeated, usually periodically at the JC Rabi rate 
Hit Extensive overviews are given in 

However, a different approach to the rebirth issue is 
more rewarding. Suppose one applies a modest-strength 
quantum field in each cavity, sufficient to dominate the 
JC sequence of ±1 photon exchanges with the atoms. As 
mentioned, the local atom-photon entanglement in each 
cavity has been examined in this case [!, 0] , but the delo- 
calized atom-atom entanglement has so far been ignored. 
It is potentially much more valuable in a variety of quan- 
tum information storage and/or transport applications. 



II. EVOLUTION EQUATIONS 

Quantum field irradiation with separate coherent 
states is known to induce very complicated AB two-party 
dynamics [2(| . Here we analyze this simplest coher- 
ent but quantum case. With a calculational trick we 
obtain a formula that suggests the surprising possibility 
of on/off switching of non-local entanglement over times 
of operation that are much longer than the Rabi cycle 
time. A direct consequence will be the ability to "hide" 
entanglement for a substantial time and recover it deter- 
ministically, and if not used hide it again, repeatedly. 

The description of CQED evolution is via the Jaynes- 
Cummings [2] interaction, which is governed by the fa- 
miliar Hamiltonian (with h — 1): 

fftot = + 5( at(7 - + a + a ) + wata 

+Y^+9(b^+^b)+wb\ (1) 

where ujq is the transition frequency between the two lev- 
els of the atoms, g is the constant of coupling between the 
atoms and the fields and u> is the angular frequency of the 
single-mode field. The usual Pauli matrices describe the 
atoms, while a', a and b\ b are the raising and lowering 
operators for the fields in the two single-mode cavities. 
A case in which the frequencies and detunings are not 
equal in the cavities has been examined [2l| . 

The JC eigenstates, as superpositions of the bare atom 
and cavity product states \g;ri) and |e;n — 1) are well 
known In order to calculate the time evolution of 
a joint atom-atom state we need to first calculate the 
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time evolution of the states of the individual sites, and 
for either site A or site B we have: 
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|e;n) = cos(<7<V n + 1) |e; n) 

— i sin(gty/n + l)|g; n + 1) 
= cos(gtVn)|.g;n) 

— i sin(gt\/n)\e; n — 1). 



(2) 
(3) 



III. TWO-QUBIT THEORY WITH COHERENT 
STATE FIELDS 



We assume AB quantum information has been stored 
(qubit entanglement has been arranged) prior to t = 0, 
for example in the pure Bell State 



|^ B (0)) = (|e<?) + \ge))/V2. 



(4) 



The coherent state characterized by n — \a\ 2 is given by 
l«) = Y^n^Mn), where A n = e~W 2 / 2 a n / V^.. Our 
initial state for the whole system is therefore 



|¥tot(Q)> = |*ab(0)) ® |a) 



(5) 



Using these results, the time evolution of the initial 
state of the system is found to be given by the double 
sum, 
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|* 4ot (t)) = — ^ ^ 



(0) 



where A^ mra is given by the formidable expression 

K m n = ~iC n+ iS m \e, e;n, m - 1) + C n +\S m \e,g; n,m) 

- S n+ iS n \g 1 e:n+ l,m) - i5 n+ iC n |p, g; n + l,m) 

- iS n C m+ i\e,e;n-l,m+l) 

- SnS m+ i\e,g;n — l,m + l) 
+ C„C m +i\g,e;n,m + l) 

- iC n S m +i\g,g;n,m + 1). (7) 

with the abbreviations C n = cos^i-^/n) and «SVi = 
sin^i-y/n). 

By tracing the photon states from \^tot(t))(^tcyt{t)\ we 
obtain the 4x4 reduced density matrix pab for the two 
atom qubits, whose entanglement we will follow. Be- 
cause of the infinite range of photon numbers in a coher- 
ent state, this density matrix is drastically different from 
the mostly-zero A-matrix [22j found in almost all prior 
discussions of rebirths, becoming a matrix with no zero 
elements at all. That is we have: 



PAB 



a 



w* 



w 

b z 

z* c 

d 



a x x x 

x b z x 

x z* c x 

x x x d 



(8) 



However, by adopting a trick described below, whose va- 
lidity has to be checked numerically, the elements of pab 



marked x can all be set to zero. This doesn't eliminate 
the doubly infinite sums, but it provides a simplification 
sufficient to lead to a relatively compact final formula, as 
follows. 

We adopt Wootters' concurrence C [H| (1 > C > 0), 
where C = indicates separability (zero entanglement) 
and C = 1 means maximal pure state entanglement, as 
in a Bell state. The concurrence of an A state like © 
with x = everywhere is given by the simple expression 



C = 2 max{0, 



d}. 



(9) 



A different approach to measuring the entanglement, via 
untraced pure states, gives analogous information (see 
0). 

The elements z, a, d of pab are given by doubly infi- 
nite summations. For z one finds 

2 = 2 { ^ A 2 n A 2 m C n C n+ iC m C m+ i 



A n A n - 


iA m A m +iS n C n - 


-iC m S r 


+ A n A n - 


2A m A m j r 2S n S n - 


lSm+1 


— A n A n - 


iA m A m +iS n C n - 


-\S m +l 



and the series summations for a and d are; 

1 oo 

n - -I V* a 2 a 2 r 2 s* 2 



A n A n +iA m A m _iS n -\-iC n -\-iS m C m 

A 2 A2 q2(~i2 

A n A n _ i A m A m+ 1 S n C„ S m +i C m +i | (11) 



and 



d = 



1 OO 

-/ 4 2 /4 2 s* 2 r 2 

2 1 ,7^ ™ m "+ 1 ™ 

+ 

+ A n A n -\A m A m+ xS n C n S m j r iC m+ \\ . ( 
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The infinite extent of these summations reflects the fact 
that we have coupled the qubits to an open state space. 

The sums cannot be evaluated in closed form, but Stir- 
ling's formula, n\ — \/2™n" e~" , and Euler's formula for 
approximating summations by integrals, can be used for 
coherent states that are even only moderately strong, i.e., 
a > 10. The dominant contribution near n because of 
the Poisson-peaked nature of A n also justifies the approx- 
imation 



Vn+T = y/n ■ 



1 



(13) 



With these approximations, and the saddle point 
method of integration, highly simplified expressions for 
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the sums can be found [25|]. Abbreviating r = gt, we 
find, for example, 



J(r) = 



exp 



s ir jla. 



(14) 



A second integral is similar, but with exp(ir jls/n) re- 
placed by exp(2i-\/m"), and the saddle point method is 
again appropriate, although evaluation is more compli- 
cated. Helpful cancellations [25[ can be identified and 
lead to the following expression for Q(t) = \z\ — \fad~: 



Q{ f ) - 1 1 ( '-M> 



16a 4 



-r 2 /2 



cos(4ar) 



y — 



2nk 



exp 



2(r - 2nkaf 



fe=l,2,... 

x cos[4q;(t — 27rfca)] 



l+7T 2 fc 2 



(15) 



In writing this last equation we have used the fact that 
around t = 2nka only the term with the correspond- 
ing k gives a significant contribution to the sums. The 
contribution to r — 2irka from any other k' is propor- 
tional to exp{— Air 2 a 2 (k — fc') 2 /[l + (irk 1 ) 2 ]}, so it decays 
exponentially with the distance from k. 



IV. ANALYSIS AND SUMMARY 

There is a substantial amount of analysis behind for- 
mula (fT~5|) , and its feasibility relies on the trick that per- 
mits the step x — > in ([5]) . The trick is not hard to under- 
stand. First, the relatively narrow Poisson distribution 
of photon numbers in a coherent state suggests replacing 
I a) Cg) I a) by the Fock state |n) <g> |n), since we expect all 
main effects to be concentrated in the near neighborhood 
of the Poisson peak at n ~ n — \a\ 2 . Second, we assumed 
that the coherent states were close enough in amplitude 
to use the same parameter a for both. 

The combined effect of the local fields induces growth 
of the elements a and d, which are the ones not already 
present in the density matrix of the original entangled 
state (UJ) . Their growth and any decline of z will cause en- 
tanglement to decrease. Inspection of formula (fT5j) shows 
that it contains repeated zeros for Q(t) — \z\ — Vad, 
meaning repeated deaths and rebirths of entanglement. 
But between clusters of death and rebirth (TT5|) predicts 
substantial intervals of time when Q remains zero to a 
very good approximation. This is demonstrated in Fig . 
[TJ where we display concurrence evolution plots (see [20|) 
from both the trick formula and the corresponding nu- 
merical solution to the evolution, without using the trick. 
Note that concurrence is plotted over a time interval 
much longer than a Rabi period: gAr = it. That is, 
we have entered the revival time regime explored earlier 
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FIG. 1: The analytical and numerical results for entangle- 
ment. The former are for the X-type p while the numerical 
ones are for the original p. The shape, location and strength 
of the revivals are predicted well by analytical formula (|15[) . 



[H 9 , but now with results for non-local two-cavity en- 
tanglement. 

Clearly, both curves in Fig. [T] show revivals, i.e., re- 
peated eruptions of finite concurrence, separated by rela- 
tively long intervals of no appreciable entanglement. This 
is the specifically quantum effect of the coherent states, 
acting here as reservoirs in the joint state evolution. As 
with all quantum revival effects, this is due to the granu- 
larity of the coherent state, to its basis in discrete photon- 
number states (fractional photons do not exist). This as- 
pect has been demonstrated experimentally several times 
(most recently, see |26j,|27|). The two plots are in remark- 
ably good agreement as to position and height of revivals, 
despite the crude trick played to get the first plot. 

The location of the rebirths is controlled by a in the 
usual way of revivals, and the kth peak height Hk(r) 
decreases with k as given by the formula 



2H k (r) - 2/7rfc - 1 + exp(-r 2 /16a 4 



(16) 



In Fig. [21 which shows an expanded snapshot of the 
curves near to the revival at r = 207r, the analytical 
formula is seen to predict a micro-structure of repeated 
ESD events on the Rabi scale. This is not present in 
the numerically exact plot, which shows no intra-revival 
death-birth events, and this is true even in the revival 
immediately following t = 0. 

Figures [T] and [5] demonstrate what we have claimed, 
that the quantum nature of the reservoirs, their photonic 
granularity, can create a new and much longer time scale 
for management of non-local entanglement. Addition- 
ally, on that longer time scale, which is many times the 
Rabi cycle period, a remarkably complete control of the 
joint non-local atom-atom state can be achieved. This 
complete control amounts to an on/off switch for entan- 
glement. Between the revivals entanglement can be con- 
sidered entirely hidden, but almost fully preserved for 
recovery at the next revival. 
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FIG. 2: A detailed plot for a = 10 of the results around 
t = 20n/g. Analytical results are for the Y-type p while the 
numerical ones are for the original p. 

To test that these conclusions are valid generally, and 
not only for special values of coherent field strength, we 
show a test in Fig. [3] for two smaller values of mean 
photon number: n = 25 and 36. In those graphs the 
recovery times are different, and even for those nearer- 
vacuum fields they still illustrate event-timing control, 
and one sees the same compact time-zone of non-zero 
concurrence through almost all of the revival episodes. 
The modulations on top of the curves can be analyzed, 
but the key is that they are small and are smaller for 



higher n, as in Fig. [2j In sum, both timing and relative 
smoothness in C are controllable features. 
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FIG. 3: Similar to the lower curve in Fig. [2] above, except 
that here a — 5 and 6, rather than 10. 



V. ACKNOWLEDGEMENTS 

We are pleased to thank Prof. Ting Yu for consultation 
and collaboration in the early stages of this study. Partial 
financial support was provided by grants from DARPA 
HR0011-09-1-0008 and ARO W911NF-09-1-0385. 



[1] T. Yu and J.H. Eberly, Phys. Rev. Lett. 97, 140403 
(2006). 

[2] E.T. Jaynes and F.W. Cummings, Proc. IEEE 51, 89 
(1963). 

[3] J. Gea-Banacloche, Phys. Rev. Lett. 65, 3385 (1990). 
[4] S.J.D. Phoenix and PL. Knight, Phys. Rev. A 44, 6023 
(1991). 

[5] J.H. Eberly, N.B. Narozhny and J.J. Sanchez- 
Mondragon, Phys. Rev. Lett. 44, 1323 (1980). 

[6] N.B. Narozhny, J.J. Sanchez-Mondragon and J.H. 
Eberly, Phys. Rev. A 23, 236 (1981). 

[7] H.-I. Yoo and J.H. Eberly, Phys. Rep. 118, 239 (1985). 

[8] R. W. Rendell and A. K. Rajagopal, Phys. Rev. A 67, 
062110 (2003). 

[9] T. Yu and J.H. Eberly, Phys. Rev. Lett. 93, 140404 
(2004). 

[10] M. Ikram, F.L. Li and M.S. Zubairy, Phys. Rev. A 75, 
062336 (2007). 

[11] M. Yonag, T.Yu and J.H. Eberly, J. Phys. B 39, S621 
(2006). 

[12] K. Zyczkowski, P. Horodecki, M. Horodecki, and R. 

Horodecki, Phys. Rev. A 65, 012101 (2001). 
[13] S. Daffer, K. Wodkiewicz, and J. K. Mclver, Phys. Rev. 

A 67, 062312 (2003). 
[14] L. Diosi, in Irreversible Quantum Dynamics, edited by F. 

Benatti and R. Floreanini (Springer, New York, 2003), 

pp. 157-163. 



[15] P.J. Dodd and J.J. Halliwell, Phys. Rev. A 69, 052105 
(2004). 

[16] A number of treatments include S. Chan, M. D. Reid 
and Z. Ficek, J. Phys. B 42, 065507 (2009), I. Sainz 
and G. Bjork, Phys. Rev. A 76, 042313 (2007), N.S. 
Williams and A.N. Jordan, Phys. Rev. A 78, 062322 

(2008) , C.E. Lopez, G. Romero, F. Lastra, E. Solano 
and J.C. Retamal, Phys. Rev. Lett. 101, 080503 (2008), 
M. Yonag, T.Yu and J.H. Eberly, J. Phys. B 40, S45 
(2007), and X.F. Qian and J.H. Eberly, in preparation. 

[17] T. Yu and J.H. Eberly, Science 323, 598 (2009). 

[18] S. Das and G.S. Agarwal. larXiv:0905.339"9l . 

[19] M. Yonac and J.H. Eberly, Opt. Lett. 33, 270 (2008). 

[20] M. Yonac and J.H. Eberly, submitted to Phys. Rev. A 

(2009) . 

[21] See S. Chan, M. D. Reid and Z. Ficek, J. Phys. B 42, 
065507 (2009) and I. Sainz and G. Bjork, Phys. Rev. A 
76, 042313 (2007). 

[22] T. Yu and J.H. Eberly, Quant. Inf. Comp. 7, 459 (2007). 

[23] W.K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[24] X.F. Qian and J.H. Eberly, in preparation. 

[25] M. Yonag, Ph.D. Thesis, Department of Physics and As- 
tronomy, University of Rochester (2009). 

[26] D. M. Meekhof et ai, Phys. Rev. Lett. 76, 1796 (1996). 

[27] M. Brune et ai, Phys. Rev. Lett. 76, 1800 (1996). 



